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ABSTRACT 



Aims. Relations connecting gamma ray burst quantities can be used to constrain cosmographic parameters of the Hubble law at 
medium-high redshifts. 

Methods. We consider a sample of 27 gamma ray bursts to construct the luminosity distance to redshift relation and derive the values 
of the parameters q , j , and s . The analysis is compared with other methods in the literature. 

Results, gamma gay bursts, if calibrated by SNela, seem reliable as distance indicators and give cosmographic parameters in agree- 
ment with the ACDM model. 
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1. Introduction 

A class of very accurate standard candles, the supernovae la 
(SNela), has been highly developed in the last two decades 
(Branch & Tammann 1992I) : however, these objects are hardly 
detectable at redshifts higher than z = 1.7, so the study of more 
distant regions of the Universe leads to the needing to implement 
more powerful standard candles. The problem becomes particu- 
larly crucial at intermediate redshift, z = 6 - 7, where, up to now, 
not very well - defined distance indicators are available. 

In the last years, several efforts have been made in or- 
der to implement gamma ray bursts (GRBs), the most power- 
ful explosions in the Universe, as standard candles, and sev- 
eral interesting results have recently been achieved (e.g. Amati 
et al 2008; Basilakos & Perivolaropoulos 2008 and references 
therein). Considering the standard model of such objects, the 
GRB phenomenon should originate from the black hole forma- 
tion and reach huge amounts of energy (up to 10 54 <?rg). These 
events are observed at considerable distances, so there are sev- 
eral efforts to frame them into the standard of cosmological dis- 
tance ladder. 

In the literature, several more - detailed models give account 
for the GRB formation, e.g. Meszaros 2006; Ruffini et al 2008, 
but, up to now, none of them is intrinsically capable of connect- 
ing all the observable quantities. For this reason, GRBs cannot 
be used as standard candles. Despite this shortcoming, there are 
several observational correlations among the photometric and 
spectral properties of G RBs. These fea tures allow use of GRBs 
as distance indicators (ISchaefer 20071) . even when they cannot 
be fully "enrolled" in the class of standard candles. In particu- 
lar, it is possible to connect the peak energy of GRBs, E p , with 
the isotropic energy released in the burst, Ej SO , and with the rest 
frame jet break - time of the afterglow optical light curve, mea- 
sured in days, %, dLiang & Zhang 20051) : 



E p (l+z) t b 

log E; m — a + b\ log ¥ bj log 

s uo B 3QQkeV 2 B (l+ Z )lday 



(1) 



where a and b,, with i = 1,2, are calibration constants. 

An other interesting resul t is the relation given by Ghirlanda 
et al. dGhirlanda et al. 2 004). It connects the peak energy E p 
with the collimation-corrected energy, or the energy release of 
a GRB jet, E 7 , where 



(1 - cos Jet )E iso , 



with 8 j et the jet opening angle, given by dSari et al. 19 99): 



"jet 



0.163 



\l+z) \E is „ 52i 



(2) 



(3) 



where E; JO) 52 = E,j O /10 ergs, «o is the circumburst particle 
density in 1 cm 3 , and r\ y the radiative efficiency. The Ghirlanda 
et al. relation is 



log E y — a + b log 



300/teV 



(4) 



where a and b are two calibration constants. 

These two relations are used the most in constraining cos- 
mology due to their relatively small scatter, interestingly very 
tight in the Ghirlanda et al. one, and the sufficient number of 
data p oints available. 

In dSchaefer 20 07). an example of the discrepancy between 
data and theoretical curves is shown for these two relations. 
It is worth noticing that the calibration of the above relations 
is necessary to avoid the circularity problem: all the relations 
need to be calibrated for every set of cosmological parame- 
ters. Indeed, all GRB distances, obtained only in a photometric 
way, are strictly dependent on the cosmological parameters since 
there is no low-redshift set of GRBs to achieve a cosmology- 
independent calibration. 

Recently, Liang et al. (Liang et al 2008) present a calibration 
method (Liang thereafter) for several GRB relations, included 
the above relations ([TJ and ©, in a cosmology-independent way 
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Table 1. Parameter values obtained by (Lian g~et al 20081) 



Relation 


a 




b 




Ey — Ep 

Ejso ~ Ep — tfj 


52.26 ± 
52.83 ± 


0.09 
0.10 


1.69 ± 
2.28 ± 
-1.07 ± 


0.11 
0.30 
0.21 



using the SNela. In fact, the SNela are very accurate standard 
candles, but their range is limited up to z ~ 1 .7; hence, assuming 
that relations (HJ and d3]i work at any z and that, at the same 
redshift, GRBs and SNela have the same luminosity distance, it 
becomes possible, in principle, to calibrate GRB relations at low 
redshifts. The calibration parameters are shown in Table Q] For 
the Ei SO — E p — tb relation, the b- values in the first line is b\ and 
in the second line is b^- 

When our working-relations are calibrated with the Liang 
method, we can compute the luminosity distance di from the 
well-known relation between d/ and the energy-flux ratio of 
the distance indicators in consid eration. Afte rwards, we can use 
a formulation given by Visser (Visser 2004), where the lumi- 
nosity distance di is related to the cosmographic parameters 
(I Weinberg 19721) by means of a Taylor series expansion for the 
same d\. Such an analysis works very well at low and interme- 
diate redshifts, since very good classes of standard candles are 
available there. Besides, it is useful to constrain alternative theo- 
ries of gravity, as shown in Capozziello et al. 2008. Since we are 
calibrating GRBs by SNela (in the SNela redshift range, the d[ 
Taylor series analysis works very well), the method could also 
be extended to the next step (intermediate-high redshifts) where 
GRBs are expected to be suitable distance indicators. This work- 
ing hypothesis could be useful in order to link low and high red- 
shift ranges and then fully probe d/. However, it is clear that such 
a Taylor expansion, derived for low redshifts, can be problematic 
for fitting GRBs at high redshifts. Here, we consider it a viable 
methodological approach to link GRBs to SNela. 

T he aim of this wo rk is to achieve the cosmographic param- 
eters (I Weinberg 19721) using the above GRB relations and then 
to test the cosmological density parameters in a ACDM model. 
The only assumption that we make here is that the Universe is 
described by a Friedmann-Robertson- Walker geometry and the 
scale factor of the universe a(t) can be expanded in a Taylor se- 
ries (Sect.2). In Sect. 3, after considering a sample of 27 GRBs, 
we use a best-fit analysis to derive the cosmographic parame- 
ters discussed in the previous section, adopting the so - called 
Chevallier, Polarsky, Linder parameterization for the equation 
of state (EoS). Discussion and conclusions are given in Sect.4. 



To build a distance-redshift diagram, one has to calculate the 
luminosity distance for each GRB in a given sample. In our case 
the luminosity distance is 



di 



bolo 



(5) 



where S' boh - S boiolQ + z) is the bolometric fluence of gamma 
rays in the burst, corrected with respect to the rest frame. The 
definition of E, so is different for each relation used, therefore for 
the luminosity distance, we have 



di — 



1U \ 300keV j \{\+z)\day) 



4ttSI 



bolo 



1/2 



(6) 



adopting the Liang-Zhang relation, with a, b\, and b^ given in 
the Table Q] and 



, n .(1 + z)a 2 ^[E p (\ +z)/100keV] 2fe / 3 
di = 7.575 — r-T r-7 Mpc, 



(7) 



for the Ghirlan da et al. relation, with a and b given in [TJ 
dXu et al. 20051) . Note that the former gives d\ in centimeters, 
therefore it divides the result for the value of 1 parsec in cm, 
while the latter gives d[ directly in Mpc. 

The luminosity d istance can be connected to the Hubble se- 
ries dWeinberg 1972b . Expanding the Hubble law up to the fourth 
order in redshift and considering the related luminosity distance, 
we get (I Visser 2 004) 



d,(z) = d H z\ 1 + i [1 - qo] Z - - 



1 



qo ■ 



-j k d Tj 

3q 2 + j + 



Z 2 



+ — [2 - 2q - I5ql - I5q\ + 5j (l + 2q ) + s 



]z 3 +<5(z 4 )X8) 



2kd 2 H (\ + 3qo) 



where dn = c/Hq is the Hubble radius and where the cosmo- 
graphic parameters are defined as 



H(t) 



1 da 

+ . 

a dt 



q(t) = - 



1 d 2 a 
a dt 2 



1 da 
a dt 



(9) 



(10) 



2. Cosmography 

The calibration we want to achieve should be cosmologically 
model-independent; hence, applying the above relations to a 
GRB sample in a given z-range, we want to derive the related 
cosmography. In parti cular, we wan t to obtain deceleration, jerk, 
and snap parameters dVisser 20041) and compare them with the 
current values deduced by other methods and observations (see, 
for example, Basilakos & Perivolaropoulos 2008; Capozziello et 
al 2008 and references therein). 

Being only related to the derivatives of the scale factor allows 
to fit the cosmographic parameters versus the distance-redshift 
relation without any a priori assumption on the underlying cos- 
mological model but, this fails at very high redshifts where the 
Taylor expansion does not work yet. 



m = 


1 

H — 

a 


d^a 
IF 


1 

a 


da 
dt 


-3 


s(t) = 


1 

+ - 

a 


d A a 
~a¥ 


1 

a 


da 
dt 


-4 



They are usually referred to as the Hubble, deceleration, jerk, 
and snap parameters, respectively. Their present values, which 
we denote with a subscript 0, may be used to characterize the 
evolutionary status of the Universe. For instance, q^ < de- 
notes an accelerated expansion, while jo allows us to distinguish 
among different accelerating models; a positive value of jo in- 
dicates that, in the past, the acceleration reversed its sign. In 
this paper, according to the WMAP observations, we assume 
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the value of the Hu bble constant Ho — 70 + 2 km/sec/Mpc 
dKomatsu et al 20081) . 

The cosmographic parameters can be expressed in terms of 
the dark energy density and the EoS. Following the p rescriptions 
of the Dark Energy Task Force, (lAlbrecht et al 20061) . we use the 
Chevallier-Polarski-Linder parame terization (CPL) for the EoS 
(Chevallie fet al 20011; Hinder 20031) and assume a spatially fiat 
Universe filled with dust matter and dark energy. The dimension- 
less Hubble parameter E(z) = H/Hq reads as 



E 2 ( Z ) = n M (l + z) 3 + Qx(l + z) 



3(l+wo+w a ) 



e 



(13) 



with Q.x = \ - Q.m and wo and w fl the CPL parameterization 
for the EoS (see Chevallier et al 2001; Linder 2003; Capozziello 
et al 2008 for details). We can have = Qa, with A the cos- 
mological constant. Such a relation can be used to evaluate the 
cosmographic parameters, obtaining 



1 3 

<7o = - + -(1 - ^m)w ■ 



70 = 1 + ^ ~ ^m) [3w (1 + w ) + w a ] , 



(14) 
(15) 



and substitute the equation for d\, we obtain a logarithmic ver- 
sion of the Hubble series: 

ln[di/(zMpc)] = \n(d H /Mpc) - ^[-1 + q Q ]z 

1 kd 2 
+ ^rt-3 + IO40 + 9q 2 - 4(jo + 1 + -f)]z 2 



+ 2^;[4?o(;o + 1 + kd 2 H a 2 ) + 5-9q - 16q 2 Q - I0q 3 



+7o(7 + 4« 7 „) + io ]z 3 +0(z 4 ). 



(21) 



This logarithmic version shows the advantage when there is 
no need to transform the uncertainties on the distance modu- 
lus. With these considerations in mind, we perform a polynomial 
least-squares fit for each relation of the data assuming Taylor se- 
ries polynomials, both in distance and in logarithmic distance. 
We stop at order n — 3 both for the polynomial fit and for the 
logarithmic fit. In the latter case, we obtain an estimate of the 
snap parameter. Note that we are using least squares since, in 
absence of any better data-fitting procedure, this is the standard 
procedure when assuming Gaussian distributed uncertainties. 

The truncated polynomial used in the fits has the form 



7 33 

SO = - -^-(1 - &M)w a 

9 

- -(1 - Q M ) [9 + (7 - Q M )w fl ] w 

- ^(l-Q M )(16-3Q M )w 2 

- ^(l-Q M )(3-n M )w^. 



(16) 



For a ACDM-universe, where (wq, w fl ) = (-1,0), it becomes 



qo = -1 + -£2 M ; 

7o = i; 



*o - 1 _ -&m > 



(17) 
(18) 
(19) 



that are the quantities which we are going to fit using a given 
GRB sample. 



3. GRB data fitting 

Let us take a GRB sample into account that satisfies the above re- 
lations. Unfortunately only 27 GRBs have observed jet breaks in 
the Schaefer sample (ISchaefer 20071) . The observational quan- 
tities of GRBs to take into account, are listed in Table [2] The 
luminosity distance for each of the relations is given by Eqs. © 
and (0, and then we obtain a data distribution in the luminosity 
distance-redshift diagram d\ - z. The errors on the data are only 
of a photometric nature and, in a first analysis, we can exclude 
errors o n the redshift. Fo r each GRB, we assume r\ y = 0.2 and 
o- n = 0, (iFrailet al. 20011) . 

Another version of the Hubble series can be used to improve 
the data fit. If we consider the equation for the distance modulus, 



p, = 25 + 



fri(10) 



ln[d t /(lMpc)] + 25, 



d(z) 



s 

1=1 



fl/Z • 



and 

3 

]n[d(z)/{zMpc)] = J] bd 



(22) 



(23) 



i=\ 



for the logarithmic fit. In the latter case, the Hubble constant 
enters as the i = 1 component of the fit. As stated above, we use 
Hq as a constraint (a prior). 

The fits can be used to estimate the deceleration and the jerk 
parameters. The logarithmic fit is better for estimating the snap 
parameter through the values of the coefficients a,- and bi and 
their statistical uncertainties. The statistical uncertainties on qo 
are linearly related to the statistical uncertainties on the parame- 
ter b\, while the statistical uncertainties on jo and so depend non 
-linearly on qo and its statistical uncertainty. It is worth noticing 
the combin ation jp + kd 2 H la % which is a well-known degeneracy 
in Eq.(l2"n> (I Weinberg 19721) . It means that we cannot determine 
jo and Q = 1 + kd 2 H la^ separately, but we need an independent 
determination of Q. to estimate the value of the jerk parameter. 

The results of the fits are presented in Table [3] and all of 
them include the error on the data. For the calculation of the un- 
certaintie s on d\, we hav e followed the procedure discussed in 
Xu et al. dXu et al. 20051) . For example, the fractional uncertain- 
ties on d\ in the Ghirla nda et al. relation , without the small angle 
approximation for 0j et (ISari et al. 1999b . are given by 



°~d L 
d L 



<rs b , 



>bolo 



1 



4(1 



3cr,,, 



I , °"i , E P 

+ {b - ln Too 



«o 



C 9 ) 2 

2 



i 

+ - - 



Co 



4(1- VQ) 2 



(24) 



where Cg = [6 sin 0/8-8 cos 8] 2 . This shows that the uncertain- 
*^0) ties on di in the Ghirlanda et al. relation are very high, due to the 
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Table 2. GRBs Data Sample 



GRB 


z 


E p (keV) 


Sboio (erg 


cm l ) 


t jet (days) 


6 je , (deg.) 


rig (cm j ) 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


970508 


0.84 


389 


±40 


8.09E-6 ± 


8.1E-7 


25 ±5 


23 ±3 


3.0 ± 


2.4 


970828 


0.96 


298 


± 30 


1.23E-4± 


1.2E-5 


2.2 ± 0.4 


5.91 ± 0.79 


3.0 ± 


2.4 


980703 


0.97 


254 


± 25 


2.83E-5 ± 


2.9E-6 


3.4 ± 0.5 


1 1.02 ± 0.8 


28.0: 


t 10 


990123 


1.61 


604 


±60 


3.11E-4± 


3.1E-5 


2.04 ± 0.46 


3.98 ± 0.57 


3.0 ± 


2.4 


990510 


1.62 


126 


± 10 


2.85E-5 ± 


2.9E-6 


1.6 ±0.2 


3.74 ± 0.28 


0.29 ± 


0.14 


990705 


0.84 


189 


± 15 


1.34E-4± 


1.5E-5 


1 ± 0.2 


4.78 ± 0.66 


3.0 ± 


2.4 


990712 


0.43 


65 : 


t 10 


1.19E-5 ± 


6.2E-7 


1.6 ± 0.2 


9.47 ± 1.2 


3.0 ± 


2.4 


991216 


1.02 


318 


± 30 


2.48E-4 ± 


2.5E-5 


1.2 ± 0.4 


4.44 ± 0.7 


4.7 ± 


2.8 


010222 


1.48 


309 


± 12 


2.45E-4 ± 


9.1E-6 


0.93 ±0.1 


3.03 ± 0.14 


3.0 ± 


2.4 


011211 


2.14 


59 


± 8 


9.20E-6 ± 


9.5E-7 


1.56 ±0.16 


5.38 ± 0.66 


3.0 ± 


2.4 


020124 


3.20 


87 : 


t 18 


1.14E-5 ± 


1.1E-6 


3 ±0.4 


5.07 ± 0.64 


3.0 ± 


2.4 


020405 


0.70 


364 


±90 


1.10E-4± 


2.1E-6 


1.67 ± 0.52 


6.27 ± 1.03 


3.0 ± 


2.4 


020813 


1.25 


142 


± 14 


1.59E-4± 


2.9E-6 


0.43 ± 0.06 


2.8 ± 0.36 


3.0 ± 


2.4 


021004 


2.32 


80: 


t 53 


3.61E-6 ± 


8.6E-7 


4.74 ± 0.5 


8.47 ± 1.06 


30.0 ± 


27.0 


030226 


1.98 


97 : 


t 27 


8.33E-6 ± 


9.8E-7 


1.04 ± 0.12 


4.71 ±0.58 


3.0 ± 


2.4 


030328 


1.52 


126 


± 14 


6.14E-5 ± 


2.4E-6 


0.8 ±0.1 


3.58 ± 0.45 


3.0 ± 


2.4 


030329 


0.17 


67.9 


± 2.3 


2.31E-4± 


2.0E-6 


0.5 ±0.1 


5.69 ± 0.5 


1.0 ± 


0.11 


030429 


2.66 


35: 


t 12 


1.13E-6± 


1.9E-7 


1.77 ± 1.0 


6.3 ± 1.52 


3.0 ± 


2.4 


041006 


0.71 


63: 


t 12 


1.75E-5 ± 


1.8E-6 


0.16 ±0.04 


2.79 ± 0.41 


3.0 ± 


2.4 


050318 


1.44 


47 : 


t 15 


3.46E-6 ± 


3.5E-7 


0.21 ± 0.07 


3.65 ± 0.5 


3.0 ± 


2.4 


050505 


4.27 


70: 


t 23 


6.20E-6 ± 


8.5E-7 


0.21 ± 0.04 


3.0 ±0.8 


3.0 ± 


2.4 


050525 


0.61 


81.2 


± 1.4 


2.59E-5 ± 


1.3E-6 


0.28 ±0.12 


4.04 ± 0.8 


3.0 ± 


2.4 


050904 


6.29 


436: 


t 200 


2.0E-5 ± 


2E-6 


2.6 ± 1 


8± 1 


3.0 ± 


2.4 


051022 


0.80 


510 


± 22 


3.40E-4 ± 


1.2E-5 


2.9 ± 0.2 


4.4 ±0.1 


3.0 ± 


2.4 


060124 


2.30 


237 


±76 


3.37E-5 ± 


3.4E-6 


1.2 ± 


3.72 ±0.15 


3.0 ± 


2.4 


060210 


3.91 


149 


± 35 


1.94E-5 ± 


1.2E-6 


0.33 ± 0.08 


1.9 ± 0.17 


3.0 ± 


2.4 


060526 


3.21 


25 


±5 


1.17E-6± 


1.7E-7 


1.27 ± 0.35 


4.7 ± 1 


3.0 ± 


2.4 



References: (Jimenez et al. 2001); (Metzger et al. 1997); (Djorgovski et al. 1999); (Kulkarni et al. 1999); (Israel et al. 1999); 

dBjornsson et al. 20011) : dLi et all) 



Table 3. Results of the fits. LZ is for Liang-Zhang relation, GGL 
for the Ghirlandaet al. one 



Fit 


qo 




jo 


+ n 




•so 


d,(z) LZ 


-0.94 ± 


0.30 


2.71 


± 1.1 






d,(z) GGL 


-0.39 ± 


0.11 


2.52 


± 1.33 






ln[d,/z] LZ 


-0.68 ± 


0.30 


0.021 


± 1.07 


3.39 


± 17.13 


ln[d,/z] GGL 


-0.78 ± 


0.20 


0.62 


±0.86 


8.32 


± 12.16 



dependence on several parameters. For this reason, in FigQ] the 
prediction bounds are plotted at a 68% confidence level, instead 
of 95%, as in LZ-relation. 

As said above, only statistical uncertainties have been con- 
sidered and other kinds of errors (systematics of cosmologi- 
cal inference, modelling errors and "historical" biases, Visser 
2007b) have been neglected. If we do not assume Ho as a con- 
straint, the analysis gives Ho = 56 km/s/Mpc, which means that 
the data sample needs to be improved with further GRBs to give 
more reliable results. 

Another step would be to test the goodness of the next fit 
statistics using the MATLAB package. In particular, we used the 
R-square method: A value closer to 1 indicates a better fit. In 
Table |4] the results of R-square are shown and the plots of the 
residuals of the fits are shown in Figs. [T] [2] For the logarithmic 
fit, the bad value of the R-square is caused by the logarithm of the 
Hubble series, which spreads a lot of the data on the /«(ii;)-axis. 
The values <K 1 for the logarithmic fits are due to the discrepancy 
of the data. 

In summary, the results are in quite good agreement with 
the ACDM model, giving a Universe model that accelerates in 
the present epoch and that has undergone a decelerated phase in 



Table 4. Goodness of the fits with the R-square. 



Fit 


R-square 


d,(z) LZ 


0.9909 


d,(z) GGL 


0.9977 


ln[d,/zMpc] LZ 


0.4005 


ln[d,/zMpc] GGL 


0.2929 



the past. The signature of this past phase is related to the sign 
change of the parameter qo and the positive value of the jerk 
parameter, unless a positive value of the spatial curvature con- 
stant k is considered. However this occurrence is excluded by the 
last observational res ults, which confirm a spatially flat Universe 
dKomatsu et al 20081) . 

3.1. The CPL parameterization to test the ACDM model 

As said, the cosmographic parameters may also be expressed in 
terms of the dark energy density and EoS parameters. Starting 
from the Friedmann equation, we obtain the Hubble parameter: 



where p is the energy density. The continuity equation for 
each cosmolog ical component is given by the Bianchi identity 
dWeinberg 1972b : 

^ = -3//( 1 + l - ) = -3H[\ + w(z)l (26) 
P \ PI 

where p is the pressure of the component considered and w(z) = 
pip the redshift-dependent EoS for each component. The dark 
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Table 5. Cosmological density parameters 



x 10 



Data and Fit 



Fit 






d,(z) LZ 


0.04 ± 


0.03 


0.65 ± 


0.73 


di(z) GGL 


0.46 ± 


0.43 


0.54 ± 


2.82 


MdJzMpc] LZ 


0.37 ± 


0.31 


0.63 ± 


1.13 


\n[d,lzMpc] GGL 


0.28 ± 


0.30 


0.72 ± 


1.09 



energy component responsible for the obse rved acceleration 
of the universe m ust have a negative EoS, dRiess et al. 2 004; 
lAllen et al. 20041). To find w(z)de, we can adopt the CPL param- 
eterization, ( Chevallier et al 2001b iLin der 2003), where 



w(z)de = w + w a z 



1 



1 +z 



(27) 



with wo and w a two parameters that enter directly into the 
equations for the cosmographic parameters, ( 114H16I I. To test the 
ACDM model, we have assumed (wo, w fl ) = (-1,0). Conversely, 
such a case can be generalized by deducing the value of cos- 
mographic parameters (Q«, wo, w a ) from polynomial fits where 
GRB data are considered. Adopting our GRB sample, we ob- 
tained the following values 



w 



-0.53 + 0.64 



0.59 ± 0.77, 



(28) 



which directly enter into the CPL parameterization. The errors 
on the CPL parameters are directly connected with the errors 
on the cosmographic parameters, as easily seen from the system 
( fT4lfT6T ). The values of wo and w fl agree, within the errors, with 
the ACDM model, without assuming constraints a priori on the 
cosmological model. 



4. Discussion and conclusions 

Starting from some relations connecting the observable quanti- 
ties of GRBs, we have used a sample of 27 GRBs to derive the 
luminosity distance - redshift diagram of the Hubble law. The 
relations conveniently calibrated by SNela to make them inde- 
pendent of any cosmological models. 

We have taken the Hubble law into account in the Taylor se- 
ries form, assuming the luminosity distance d/ as a redshift func- 
tion whose coefficients are combinations of the cosmographic 
parameters Ho, qo, jo, and so- The aim was to evaluate such pa- 
rameters starting from the GRB data. A direct analysis of the 
fits leads to the conclusion that, in the error ra nge, the SNel a 
results can also be extended at higher redshifts (iVisser 2007bh . 
Besides, such results agree with the ACDM model according to 
Eqs. (fT7l),(fT8Tl,(fT9l). In particular, the value of the parameter qo 
that we found is in agreement with the observed Q.m (see Table 
|5]>. 

However, the sample we used is quite poor at high redshifts 
and, in some sense, this justifies the use of the method of Taylor 
series which works very well at low redshifts. In particular, at 
z > 6, we only have one GRB, GRB050904 (see Fig.Q). This 
GRB is very important in the fit results because it affects the 
trend of the fits. For this reason we need some richer sample at 
medium-high redshifts to constrain the results better. However, 
if we had richer samples at high redshifts, the Taylor series anal- 
ysis would fail to constrain cosmological models since an exact, 
and not approximated, d/(z) expression is needed in that case. 
The best constraint, however, would be an absolute relation be- 
tween several, GRB observables which would make the GRBs a 
powerful standard candle at intermediate-high redshift. 




Redshift z 



Fig. 1. Luminosity distance - redshift diagram and the residuals 
of the di(z) GGL fit. Note the discrepancy at high-z. The dotted 
lines are the bounds predicted at 68% confidence level. 




Fig. 2. Logarithmic version of the luminosity distance relation 
versus redshift and the residuals of the ln[d//zMpc] LZ fit. In this 
version of the Hubble series, the discrepancy is the same at every 
Z. The dotted lines are the bounds predicted at 95% confidence 
level. 



Considering these preliminary results, it seems that cosmog- 
raphy by GRBs could be a useful tool in constraining self- 
consistent cosmological models even if, up to now, GRBs are 
not standard candles in the proper sense. 
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